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A TECHNIQUE FOR CONSTRUCTING SPECTRAL REFLECTANCE CURVES 
FROM VIKING lANLER CAMERA MULTISPECTRAL DATA 


by 

Stephen K. Park, Friedrich 0. Huck, 
and Brent D. Martin 

SUMMARY 

This report presents and evaluates a technique for the cons cruction of 
spectral reflectance curves from multispectral data obtained with the Viking 
lander cameras. The multispectral data is limited to 6 channels in the 
wavelength range 0.4 to 1.1 jim, and several of the channels suffer from 
appreciable out-of-band response. Previous analysis showed that curve fits 
through 6 spectral data points plotted at distinct wavelengths often did not 
reveal important spectral features present in the materials and sometimes 
generated false spectral features. 

The technique proposed herein represents the estimated reflectance 
curves as a linear combination of known basis functions with coefficients 
determined to m^- s the error in the representation, and it pemits all 
channels, with and ,iout out-of-band response, to contribute equally valid 
information. The technique is evaluated for known spectral reflectance curves 
of 8 materials felt likely to be present on the Martian surface. The 
technique provides an essentially exact fit if the reflectance curve has no 
pronounced maxima and minima. Even if the curve has pronounced maxima and 
minima, the fit is good and reveals the most dominant features. Since only 6 
s^miples are available some short period features are lost. This loss is 
almost certainly due to undersampling rather than out-of-band channel response . 



INTRODUCTION 


The Viking lander cameras will image the Martian surface with 6 spectral 
channels in the wavelength range OA to 1.1 |im. It is desired to construct 
approximate spectral reflectance curves from this data to aid in identifying 
surface materials . This objective is rendered difficult not only by the small 
number of channels, but also by the appreciable out-of-band response of several 
of these channels . 

The out-of-band response of spectral channels makes it inappropriate to 
generate data points by associating spectral reflectance values with distinct 
wavelengths . That is, it is unlikely that data points so constructed will lie 
on the true spectral reflectance curve. Consequently, any method of fitting a 
curve to these data points is quite unlikely to adequately approximate the 
true spectral reflectance curve. A preliminary evaluation (ref. l) confirmed 
this suspicion. It was found that a curve fit through these 6 data points 
often did not reveal important spectral features and sometimes generated 
spectral feat'ures which were false. 

The purpose of this paper is to present and evaluate a more appropriate 
technique for constructing Spectral reflectance curves. In this technique, 
the curves are represented as a linear combination of known basis functions 
with the coefficients detennined to minimize the error in the representation. 
Two sets of basis functions are considered which generate, respectively, a J^th 
degree polynomial and a cubic spline (with 6 knots). The technique is eval- 
uated by using it to approximate the knovm spectral reflectance curves for 8 
materials felt likely to be present on the Martian surface. 
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INTRODUCTION 


The Viking lander cameras will ur%ge the Martian surface with 6 sjwctral 
channels in tlie wavelength range O.U to 1.1 ^m. It is desired to construct 
approximate spectral reflectance carves from this data to aid in identifying 
surface materials. This objective is rendered difficult not only by the small 
monber of channels, but also by the appreciable out-of-band response of several 
of these channels. 

The out-of-band response of spectral channels makes it inappropriate to 
generate data points by associating spectral reflectance values with distinct 
wavelengths. That is, it is unlikely that data points so constructed will lie 
on the true spectral reflectance curve. Consequently, ary method of fitting a 
curve to these data points is quite unlikely to adequately approximate the 
true spectral reflectance curve. A preliminary evaluation (ref. l) confirmed 
this s'jsplcion. It was found that a curve fit through these 6 data points 
often did not reveal Impcrtarit spectral features and sometimes generated 
spectral features which were false. 

The purpose of this paper is to present and evaluate a more appropriate 
technique for constructing spectral reflectance curves. In this technique, 
the curves are represented as a linear ccmblnation of known basis functions 
with the coefficients determined to minimize the error in the representation. 
Two sets of basis functions are considered which generate, respectively, a '^th 
degree polynomial and a cubic spline (with 6 knots). The technique is e\'al- 
uated by using it to approximate the known spectral reflectance curves for 8 
materials felt likely to be present on the Martian surface. 
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FORMULATION 


Iiet p(X) be the spectral reflectance of the unknown material that Is 
Imatjed by the camera. This reflectance can be represented as 

hj(X) 4 €(X) (1) 

where h.(X) are known basis functions, x. are (constant) coefficients to 
*1 J 

be determined, and c(X) is the error Introduced by this representation. 

As noted in the Introduction, two sets of basis functions will be 

evaluated . 

Knowledge of the reflectance P (X ) is limited to m(-6) spectral 
samples. Except for a channel dependent, multiplicative constant* these 
samples are given by m quantities, b^, where 



The transfer function Tj(X) accounts for solar irradiance and camera 
responsivity, and is given by 

Tj(X) - S(X)T^(X)T^(X)R^(X) 

where S(X) is the solar irradiance, t (X) .he atmospheric traramittance, 
T^(X) the camera optical throughput, and Rj^(X) the photosensor array 
responsivity for the I'th chemnel. Figure 1 presents plots of Tj^(X) 

* Each of these channel dependent constants are determined by a cali- 
bration using as reference a test chart (ref. l) containing surfaces with 
known spectral reflectance curves. 
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for the 6 narrowband channels of the Viking lander camera. 

In practice, must be estimated from digital data that is contaminated 

by electronic and quantization noise as described in referenje 2. However, 
the resulting uncertainties in the actual value of b^ are not considered 
here. 

Substituting equation (l) into (2) yields m equations 


b 


i 



* "i 


i ■ 1,2, ... , ffl 


( 5 ) 


where 


‘ij ‘jf* 


and 


*i 



€(X)T^(X )dX 


(M 


{•?) 


The quantities a^^ can be calculated 
hj (X ) has been made. The error terms 
between the unknown reflectance p(X) 


once a choice of basis functions 

provide a measure of the d ‘ ■'i • * -?nce 
r 

and the estimate ]P x. h.(X). 

>1 ^ ^ 


If e^ ■ 0 for i • 1,2,..., m, then, as measured by the 6 spectral channels, 

the reflectance and the estimate would be indistinguishable (i.e., the 

reflectance and the estimate would generate identical samples b^, b^, ..., 

b ). However, this does not mean that the estimate and the unknown reflectance 
m 

are necessarily equal for all wavelengths (i.e., that €(X)»0). 

In general, the unknown coefficients Xj should be chosen to minimize 


U 



the luin of squared errors; that is. 


♦ e. 


e 


D 


min 


or equivalently the root -me an- square (ms) ermr 



-» min 


where from equation (5) 


e 


i 


“ b 


i 



( 6 ) 


This minimization problem is the well-known (linear) least squares problem, 

and it can be solved on a digital computer by standard algorithms. 

Two special, cases might help to illustrate the technique that is 

presented herein. Let the number of spectral samples be equal to the number 

of coefficients (i.e., m ■ n ■ 6). If the basis functions ^j(X) are chosen 

so that the square matrix A with elements a. . is non-singuxar (i.e., the 

* 0 

inverse, A” , exists), then ej^ ■ 0 for i ■ 1,2,..., m and the column 
vector X with components Xj^,Xg, . . .,x^ is uniquely detemined by x ■ A"H> 
where b is the column vector with components b^,b 2 « . . .«b^. 

That is, in this case the 6 coefficients x^,X 2 « • • are detemined by 

solvi’^g 6 (simulteuieous ) linear equations. Moreover, the reflectance curve 
euid tlie estimate are indistinguishable. 

As a second case consider sampling with spectral filters so narrow that 
it is appropriate to attach spectral reflecte...,e values to distinct wave- 



lengths. That Is, the transfer functions Tj^(X) can be approxlipated by 
impulse ‘ .'nctlons at distinct wavelengths as given by 

T^(X) - 6(X-X^), 1 . 1,2, ...,m 

and shovti in figure 2(a). Equations (2), (^i), and (6) become, respectively. 


bj - I P (X) 6 (X - X^)dX • p(X ) 

•'o 


'1 ■ 


The spectral samples bj^ and elements a^^j are the actual values of the 
surface reflectance P (X ) and basis functions h.(X) respectively at 

J 

distinct wavelengths X . The coefficients x are determined so as to force 
n “ 

the estimate Xj best fit the reflectance at the wavelengths 

X^, as is illustrated in figure 2b. The error term e^^ accounts for the 
amount by which the estimate fails to fit the reflectance at the wavelengths 


Xi. 


To summarize, the technique proposed herein for constructing spectral 


reflectance curves from m samples bj^,bg, . . .,b^ (eq.(2)) is as follows: 

(i) select a set of n basis functions h.(X) (n < m); 

iJ 

(ii) evaluate the mn quantities a^^j given by equation (U); 

(iii) determine the n coefficients Xj by either solving the linear 
least squares problem, or, if n=m emd the square matrix A with 
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eleaenti Is non-singular, by solving the linear equations (3) 

with e. - Oi 

n 

(Iv) fom the estimated spectral reflectance x. h.(X). 

>1 J J 

It should be recognized that this technique permits all channels, with and 
without out-of-band response, to contribute equally valid Information because 
it does not require spectral reflectance values to be associated with distinct 
wavelengths . 

The relative success (or failure) of the technique Is determined by the 
choice of basis functions. Certainly, some choices are more appropriate than 
others. In the next section, two sets of basis functions are discussed. 
Following that, these basi** sets are evaltiated to demonstrate how accurately 
they generate approximations to th<> knowit spectral reflectance curves of 8 
materials . 


EAT IS FUNCTIONS 

If the basis functions are chosen to be hj^(X) ■ 1, h2(X) ■ X,..., 
hn(X) ■ X*^“^ then the estimate x^X*^"^ is an (n-l)^^ degree polynomial 

In the variable X with coefficients Xj . With n=m*6 the square matrix 
A with elements 

is non-singular so that the polynomial Is uniquely determined by solving 
equations (3) with e^^ ■ 0. It is known that high degree polynomial approxi- 
mations can possess undesirable oscillations. However, the results shown 
later with n equal to 6 do not exhibit excessive oscillations (except 
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possibly at the endpoints). 

A BK>re modem and usually better method for fitting curves to data 
involves the utte of a cubic spline. A cubic spline is a curve defined on a 
sequence of subintervals. The end points of these subintervals are known as 
knots. On each subinterval the curve is a ». polynomial and at the knots 
the adjacent polynomials arc Joined in such a fashion that the resulting curve 
(the cubic spline) is continuous with continuous 1st and 2nd derivatives. 

The number of ki>ots (n) is less than or equal to the number of samples 
(m<^ ), and, for convenience, the knots (X^) are chosen to be equally spaced. 
The knots used for all the calculations discussed herein are 

Vj ■ .U25 (J-l)A J - 1,2, ...,6 


where the spacing is ^ > .12^ ^m. The basis functions are 
hj(X) - C(X-X^) j - 1,2, ...,6 

where C(X) is the function defined by (see figure 3«) 


C(X) 



+ 5A^(A-|\| ) ♦ 5A(A- |X| )2 . 5 (A_ jx] )5| |x|<d. 


(2A -|x| )^ 


A< |x| < 2 A 



On each subinterval (of width A) C(X-X ) is a cubic polynomial in X, 

d 

and at each knot C(X-\^) is continuous with continuous 1st and 2nd 
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n _ 

d«ri' iti ei. Conaequently, the estimate £ x. C(X-X ) axso has these 

1.1 J J 

properties (i.e.. It Is a cubic spline). The 6 basis functions are indicated 
in figure (5b). 

With n>tt«6 the matrix A with coefficients 


• y " C(X-Xj) T^(X)dX 


is non- singular so that the cubic spline 
determined by equations (5) with e^ • 0. 


X C(X-X ) 
J-1 ^ 


is uniquely 


m:sui.T8 

Figure 4 presents the spectral reflecta/ice curves for 0 materials that 
can be reasonably expected on the surface of Mars. Using the technique 
iscussed previously, ^th degree polynomial and cubic spline estimates were 
constructed for each material. The results are presented in figures 5 and 6 
for the polyromials and cubic splines respectively. In these two figures the 
estimate is denoted as 

<P (X) > » ^ (n - 6) 

and shown as a heavy curve; the true reflectance p(X) (from figure 4) is 
shown as a light curve, llie error between real and estimated spectral 
reflectance 


€ (X ) - p (X ) - < p (X ) > 
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li Indicated in imall inserts. 

For coenjutational purposes the 6 transfer functions Tj(X) and 8 
spectral reflectance curves p(X) were measured at 29 equally spaced wave* 
lengths X « .4, .42^, .4*>0, 1.07*v, 1.1. The resulting data arrays were 
used as inputs to determine the samples (eq.(2)) and quantities a^^ 

(eq. (4)) by n’snerlcal integration (Simpson's rule). 

For both the polynoetial and cubic spline estimates there is a tendency 
fer |<(X)| to become large rapidly as the end points (X <■ 0.4 and X ■ l.l) 
are approached. For this reason and due to people's propensity for extrapo- 
lation, estimates near the end point mast be handled with care. Thus, 
additional work is needed: 

(i) to assert a priori a reduced domain (e.g., 0.5 < ^ < 1*0) in which 
all estimates are accurate; 

(ii) or, to Impose natural boundary conditions (e.g., that the 2nd 

derivative be tern at tlie end points) in an attempt to avoid having 
to reduce the domain. 

For the 8 materials used herein all estimates are observed to be valid 
in the reduced domain .425 5^5 1*025. For those materials with relatively 
simple spectral reflectances (i.e., average Mars, Syrtls Major, pinacates, aj;d 
augite) the estimates are excellent for both the polynomial and cubic spline. 

For those materials with more complex spectral reflectances (i.e., 

I* 

llmonlte, hypersthene, and olivine) the estimates are very good. The dominaiit 
features are reproduced by both estimates; however, short period features are 
lost. This loss is almost certainly due to the small number of spectral 
channels rather than the out-of-band response of some channels. For olivine 
the polynomial and cubic spline estimates are nearly identical, while for 
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hyp«rithene and llaonlta the cubic spline estimate is superior. Specifically 
note that the cubic spline quite accurately predicts the doninant peak and 
absorption band for hypersthene. In general the results indicate a preference 
for the use of cubic splines as the basis functions. 

As mentioned previously the analysis herein ignores the fact that each 
b^ must be estimated from digital data that is contaminated by electronic 
and quaiitization noise. Also, the analysis assumes that each channel output 
can be accurately calibrated by reference to a test chart. Thus, additional 
work is necessary to determine in practice how accurately spectral reflectance 
curves can be reconstructed using as inputs samples b^ determined by the 
actual Viking lander camera channel outputs rather than samples determined 

theoretically by performing the integration indicated in equation 2. 


CONCLUDING REMARKS 

A techriique has been pr>^ sented and evaluated for constructing spectral 
reflectance curves from mu’.tispectral data obtained with the Viking lander 
cameras. The multispectral data If limited to 6 spectral channels in the 
wavelength range O.U to 1.1 urn, and some of the channels suffer from 
appieciable out-of-band response. The proposed technique presents the 
estimated reflectance curves as a linear coinbination of known basis functions 
with coefficients determined to minimize the error in the representation. An 
Important advantage of this technique is that it i>ermite all channels, with 
and without out-of-band response, to contribute equally valid Information. 

It properly accounts for the solar ir rad lance of the surface and the response 
of the camera without associating multispectral data with distinct waveien.Tths . 

The technique was evaluated for 8 spectral reflectances of materials 
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that can be reasonably expected on the surface of Mars. Estimated spectral 
reflectances for materials with relatively simple reflectances resulted In 
essentially texact fits to the actual curves. Estlnatod spectral reflectances 
for materials with relatively complex reflectances still resulted In very good 
fits; however, short period reflectance features were lost. The loss of short 
period features Is almost certainly the result of undersampling with the 
limited n\oober of channels rather than the out-of-band response of seme of 
these channels. 
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Fig. 1 - Normalized product of solar irradiance, atmospheric trans- 
mittemce, camera optical throughput, and photosensor array 
responsivity . 





Fig. 2 - Illastration of techriquo for spocial case, (a) Narrowband 
filters approxlmctp Impulse functions, (t ) Frror terms e 
occur 0 t dir^lnct waveleijgths \^. ^ 




Fig. 5 - Cubic spline, (a) Defining properties, (b) Pasls 
functions used for construction of curves. 








